Criticality without self-similarity: a 2D system with random long-range hopping 
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We consider a simple model of quantum disorder in two dimensions, characterized by a long-range 
site-to-site hopping. The system undergoes a metal-insulator transition - its eigenfunctions change 
from being extended to being localized. We demonstrate that at the point of the transition the 
eigenfunctions do not become fractal. Their density moments do not scale as a power of the system 
size. Instead, in one of the considered limits our result suggests a power of the logarithm of the 
system size. In this regard, the transition differs from a similar one in the one-dimensional version of 
the same system, as well as from the conventional Anderson transition in more than two dimensions. 
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Critical wave functions have been a subject of intense 
theoretical research during the last two decades (see 
Ref. [l| for a recent review). Recently, they have also 
attracted a lot of attention from experimentalists 0-3]. 
Starting with the innovative work by Wegner it has 
become customary to characterize the critical wave func- 
tions by their multifractal dimensions d q , which deter- 
mine the scaling of the moments I q of the wave functions 
with system size L: 

I q = L d {\^ n {v)\^)o:L- d ^- x \ (1) 

Here the averaging is performed over different disorder 
realizations, as well as over a small energy window. The 
power-law dependence of the moments with exponents 
d q =/= different from the dimensionality of the space d 
implies a self- similarity of fluctuations of the wave func- 
tion amplitudes on different spatial scales. 

Among all quantum disordered systems showing crit- 
ical behavior, two-dimensional (2D) systems occupy a 
special place. It is commonly believed that at the criti- 
cal point such systems possess conformal invariance and 
should be described by conformal field theories (CFTs), 
similar to conventional phase transitions in 2D. Although 
no direct mapping of any disordered quantum system 
onto a CFT is available at the moment, this conjecture 
of conformal invariance imposes certain constraints on 
correlation functions and critical exponents. The valid- 
ity of these predictions has been checked by numerical 
simulations for various 2D disordered systems [f| 0| , and 
no significant deviations have been found so far. One 
such prediction is a relation between the corner and the 
surface multifractalities Q: 

A c 9 (0) = ^A s g , (2) 

where the exponents A* = (d* — d){q — 1) are determined 
by the deviation of the corresponding fractal dimensions 
from the dimensionality of space d. To define one 
should take the position of the point r in Eq. |T]) to be in 
a corner of the system boundary with an internal angle 



9 (x = c), on the surface (x = s), or else in the bulk ( 
x = b). 

In the present work we study a random matrix model 
describing a long-range hopping of a particle on a 
two-dimensional lattice. Similar to its one-dimensional 
(ID) counterpart Q, the model undergoes the Ander- 
son metal- insulator transition. At the transition point, 
we compute the fractal dimensions perturbatively in the 
regime of strong criticality and find that 

d q (9) = -d q . (3) 

The dependence on 9 is inverse to that predicted by CFT 
in Eq. ([2]) . In the opposite regime of weak criticality, the 
scaling of the moments of the wave functions is found to 
be consistent not with the standard power-law ([T]), but 
with a logarithmic dependence 

I q oc L" 2 ^- 1 ) in".*'- 1 ) L. (4) 

Consequently, our model is a rather surprising example 
of a disordered 2D critical system, whose wave functions 
are not self-similar and are less localized than the stan- 
dard critical states. We believe that the origin of this 
peculiarity of our model lies in the long-range nature of 
the Hamiltonian Q. Therefore, our results might be rel- 
evant for other two-dimensional disordered critical sys- 
tems with a long-range interaction. One important ex- 
ample of such a system is the quantum Hall effect with 
unscreened Coulomb interactions. Below we present an 
exact definition of our model, derivations of the main 
equations and results of the numerical simulations sup- 
porting our findings. 

The model. The random matrix model considered 
here is defined by a random Hamiltonian matrix H r r i, 
which for simplicity can be taken in the coordinate rep- 
resentation: r and r' are two points in the system and 
H r r / is the amplitude of a quantum hop between them, 
random due to the disorder in the system. To make the 
matrix formulation meaningful, it is understood that r 
belongs to a spatial lattice. The Hamiltonian matrix can 



2 



be either real and symmetric (for a system with time- 
reversal invariance), or complex and Hermitian (for a sys- 
tem without it). In the standard random matrix theory, 
all the independent matrix elements H^, are taken to be 
independent Gaussian random variables with zero mean 
and variance e. This means that the typical hopping am- 
plitude is e for any two points. One can consider a more 
elaborate system in which the hopping amplitude decays 
with the distance: 



(5) 



We consider the case when the variance matrix depends 
only on the distance, i.e. a r ^ T > = a(|r — r'Q, and normal- 
ize the diagonal terms a(0) = 1. It is well-known that if 
a[x) decays at infinity faster than any power, all eigen- 
functions of H are localized. In the limit of a diagonal 
matrix a, this is a trivial remark. However, we are in- 
terested in the situation when the decay is a power-law 
beyond a certain bandwidth b: 



a(x) — 



1 



i + (!)*» 



(6) 



In ID, this model has been studied in depth 0, It 
is known that the system undergoes a metal-insulator 
transition at a = 1: the eigenfunctions are localized for 
a > 1 (with power-law tails) and extended for a < 1, 
irrespective of b. There is strong evidence that the eigen- 
functions are multifractal at a = 1, as in the usual 3D 
Anderson transition. Being ID, this model has an ob- 
vious advantage for calculations and hence it is widely 
used as a playground for investigating various features of 
the Anderson transition 0, [HI, El ■ 

A similar transition at a — 2 is believed to exist in 2D 
Q, where it is tempting to make a connection between 
it and a bona fide Anderson metal-insulator transition in 
2 + e dimensions. Based on the numerical evidence [l4j . 
it is certain that the eigenfunctions in the 2D power-law 
model are localized for a > 2 and extended for a < 2. 
The question remains: are the eigenfunctions fractal at 
a = 2? Since this was the case in ID, as well as in the 
true Anderson transition, the same was expected here. 
Moreover, some numerical experiments pointed in this 
direction, and even suggested the values of the multifrac- 
tal dimensions [15j. However, this numerical work was 
done only for moderate system sizes and, based on the 
results, it is not easy to judge whether or not the ob- 
tained scaling relation is real. To answer this question 
we calculate below I q in the limits of strong (b <C 1) and 
weak (6 3> 1) criticality. 

Strong criticality (b -C 1). In this regime the Hamil- 
tonian matrix H T T i is almost diagonal and we can find I q 
by virial expansion. The leading order result for q > 1/2, 
valid for any random matrix with independent Gaussian 
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FIG. 1: Dependence of the corner (6 — ~), surface and bulk 
fractal dimensions on q. Error bars are smaller than the point 
size. The solid line represents the theoretical prediction © ■ 



elements, reads [l6| 



I q {v) = 1-77 



r(«-i) 



V^T7+0(6 4 ), (7) 



where r\ = for orthogonal symmetry, 77 = for uni- 
tary symmetry. Let us assume that the system occupies 
a sector of a disk of radius L, with the point r as its apex 
and the apex angle 0. We replace the sum with an inte- 

gral: £ r , ^7 =8f Q L (l + ^ ~* %dr = simT 
and upon sending L 00, 



-1 L 
b- 



1 — rjb 



2 r(g 



5). 



r(«-i) 



lni + 0(6 4 ). 



(8) 



Comparing this result with the expected power-law scal- 
ing I q oc £-^(9-1) = L(l - (q - l)d q lnL + . . .) we can 
see that the presence of the In L term in © is consis- 
tent with the power-law ([lj , and the fractal dimension is 
given by 



d q {9) = r)t 



r(g-i) . 
r(g) 



0(6 4 ). 



(9) 



The expressions for the surface and the bulk fractal di- 
mensions can be obtained by taking 9 = n and 9 = 2n, 
respectively: 



d q (0) = ~d\ = Y ^ q . 



(10) 



The dependence on 9 is inverse to the conformal predic- 
tion P]), This result rules out the conformal invariance 
for our model at the critical point. In order to check our 
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thermore, A = l®diag(l, —1) is a matrix of the same size 
as Q, and j3 is a numeric prefactor. If the Hamiltonian 
H is real (orthogonal symmetry), f3 = 1 and the size of 
Q is 8 x 8; if if is complex (unitary symmetry), (3 — 2 
and Q is 4 x 4. The non-local kernel of the action is a 
matrix related to the Hamiltonian variance matrix: 



Ur 



(a 



ri 



From it we can construct the propagator n rr 
to 

1 

V' 



*f^r,ra — ^1*1 ,r2 



(12) 
according 

(13) 
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FIG. 2: Numerical results for the second moments of eigen- 
functions. The error bars are smaller than the point size. 
A fractal law |TJ would imply that for L — > oo the data 
points in the top panel would follow a straight line. The 
law I2 (x L~ 2 In" 2 L (a possibility suggested by Eq. (JTSJ with 
IT,,- « lnlnL) implies that the points in the bottom panel 
should follow a straight line. 



predictions Eqs. (|9I10[) we carried out numerical calcula- 
tions on square samples with lateral sizes up to L = 181 
for unitary symmetry. Eigenfunctions of the Hamiltonian 
matrix (J5J) (with the bandwidth b = 0.1) were obtained 
by standard diagonalization subroutines, and the ensem- 
ble average was performed over at least 5x 10 5 eigenstates 
with energies close to zero for each size L. 

The fractal dimensions d q were calculated from the 
scaling of the moments I q of the wave functions. In or- 
der to compute them we considered wave function ampli- 
tudes in the following regions of the LxL sample: i) 3 x 3 
squares (nine sites) at each corner for d q , ii) stretches of 
L/6 sites centered in the middle of an edge for d s q , and 
hi) an L/6 x L/6 square in the center of the sample for 

. The results are presented in Fig. [TJ The multifrac- 
tal dimensions d q are close to the theoretical predictions 
given by Eqs. (|9I10I) . 

Weak criticality (b 3> 1). In this limit the model can 
be mapped onto a non-linear er-model with a non-local 
action |9j]: 



S = -^r t Y,Ur,r>str(Q r Q r 

r,r' 



>^str(Q r A), (11) 



where Q r is an r-dependent supermatrix held, u is a tun- 
able parameter to be sent to zero in the end of the calcu- 



lation, and t 



7T 2 b^ 



1 is the coupling constant. Fur- 



where V ~ L 2 is the volume of the system. In both 
(|T2"j) and ([IB"]) the zero-modes are subtracted: J2 r ^ r , r ' = 
Y^r Hr,!-' = 0. Assuming a large translationally invariant 
system, we can work in Fourier space, where Eqs. (|12I13[) 
mean U(k) = l/a(k) - l/o(0), ft(fc) = 1/U(k) for k ^ 0, 
and n(0) = 0. For a = 2, 



a(k) 



a ikr cos <p 



1 



F" 



-rdr — . 
2tt 



To analyze a(k) at bk <C 1, we expand e lfercos ^ 
power series and integrate it term by term: 



d(k) 



a(0) - £ 



Vir 



1 



b 1 



(14) 



into a 



(15) 



The integral over r diverges and hence should be cut off 
at r ~ 1/fc. It then yields (6 4 fc 2 /4) ln(6fc), and hence the 
leading contribution to the propagator is 



1 



4 k 2 \n{bk)' 



(16) 



This method of computing the integral (|14p is quite 
straightforward, if not entirely rigorous. It can be 
checked, since this integral can also be computed ex- 
actly in terms of the Meijer G-function J.0j . Namely, 



d{k) = b T G> 



3,0 



0,i,i,0 



I ; at small arguments, this 



G-function behaves as G(y 4 ) = ir + 16(7 — 1 + log 2 + 
log y)y 2 + 0(y 4 ), from where Eq. ([TBI) again follows. 

This result is quite different from ID, where n(fc) oc 
Ifcl^ 1 , i.e. a pure power of k, which leads to multifrac- 
tality. By contrast, Eq. (ll6l) is a power-law modified by a 
logarithm, which results in the loss of fractal scaling ([T]) . 
Indeed, on the one hand, the scaling relation dTJ) implies 
that for small A„ 



L- 2 to-V(l- A q \nL 



(17) 



On the other hand, the perturbative calculation in the 
tr-model gives 



ocL- 2 ^- x )(l 



g(l-g) 
n 2 b 2 /3 



(18) 
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The two expressions are consistent if the propagator at 
coinciding points 

goes as InL for large L, in which case Ea. dT8)) yields the 
leading order expression for A q . This is the situation in 
ID. In 2D, however, Eq. (|16|) gives a double logarithmic 

2 

leading dependence TL r r m \- In In L, thus ruling out the 
possibility of the fractal scaling ([!]). Eq. with 




is the most natural law whose perturbative expansion is 
consistent with the double logarithm found above. In 
fact, one can extend the perturbative series in Eq. (|T5|> 
to the next order [Tljj , and find that it is still compatible 
with Eq. (Q}. This means that the logarithmic scaling 
hypothesis with the above value of v q holds not just in 
the leading but also in the first subleading order of per- 
turbation theory. 

In order to reconcile our result with the apparent 
power-law scaling found in (IBj , we present in Fig. [2] the 
anomalous part of I2 as a function of InL and of In InL. 
Although the system size in the simulations was at the 
limit of currently available computing power, it is still not 
very large. Hence, we do not expect the lower panel to 
represent Eq. Rather, the purpose of this panel is to 
demonstrate the impossibility of distinguishing between 
the fractal scaling law and the logarithmic one, based on 
the available data. 

In conclusion, we have demonstrated that although the 
eigenfunctions of the 2D model with power-law hopping 
at the metal-insulator transition are neither localized nor 
delocalized, they are not fractal. In the regime of strong 
criticality we have found that the scaling of the moments 
is consistent with a power-law behavior, but the relation 
between the surface and the corner fractal dimensions 
rules out a CFT description of the model. In the opposite 
regime of weak criticality, the anomalous parts of the mo- 
ments appear to scale as \aT Vq L rather than a power-law, 



making the standard notion of the wave function multi- 
fractality inapplicable. The exact form of the crossover 
between the two regimes requires further investigation. 
We believe that our results are relevant for other corre- 
lation functions as well as for different 2D models with 
long-range interactions. 
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